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ABSTRACT 
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Context. Magnetic early B-type stars are rare. Indirect indicators are needed to identify them before investing in time-intensive spectropolari- 
■ metric observations. 

" Aims. We use the strongest indirect indicator of a magnetic field in B stars, which is periodic variability of ultraviolet (UV) stellar wind lines 
I occurring symmetric about the approximate rest wavelength. Our aim is to identify probable magnetic candidates which would become targets 
^ ■ for follow-up spectropolarimetry to search for a magnetic field. 
(■ — 1 Methods. From the UV wind line variability the B1/B2V star a Lupi emerged as a new magnetic candidate star. AAT spectropolarimetric 
measurements with SEMPOL were obtained. The longitudinal component of the magnetic field integrated over the visible surface of the star 
\^ • was determined with the Least-Squares Deconvolution method. 

' Results. The UV line variations of a Lupi are similar to what is known in magnetic B stars, but no periodicity could be determined. We 
qq . detected a varying longitudinal magnetic field with amplitude of about 100 G with error bars of typically 20 G, which supports an oblique 
( ••"-) magnetic-rotator configuration. The equivalent width variations of the UV lines, the magnetic and the optical-line variations are consistent with 
£X| the photometric period of 3.02 d, which we identify with the rotation period of the star. Additional observations with ESPaDOnS attached to the 
t— I ■ CFHT confirmed this discovery, and allowed the determination of a precise magnetic period. Analysis revealed that a Lupi is a helium-strong 

t* star, with an enhanced nitrogen abundance and an underabundance of carbon, and has a chemically spotted surface. 
, , . Conclusions, a Lupi is a magnetic oblique rotator, and is a He-strong star. Like in other magnetic B stars the UV wind emission appears 
to originate close to the magnetic equatorial plane, with maximum emission occurring when a magnetic pole points towards the Earth. The 
Uh 3.01972 ± 0.00043 d magnetic rotation period is consistent with the photometric period, with maximum light corresponding to maximum 
magnetic field. 
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1. Introduction 

The star cr Lupi (HD 127381, HR 54 25, HIP 71121) ha s 
been classified as spectral type B2III by Hiltner et al. [jl969j). 
Other cla ssifications include B2V ( de Vaucouleursl 19571) and 
B2III/IV dHoukll978h . Recently. Eevenhagen & Leistetl d200rl 
concluded that B1/B2V would be more appropriate, based on 
the photosphe ric parameters eva luated in their work. 

Although Shobbrook (1978) f ound no photometric vari- 
ations, IVander Linden & Sterkenl dl9871) reported a signifi- 
cant amplitude in variations in the V band of 0.01 mag. 
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This emerged during a photometric campaign in which cr 
Lup was used as a comparis on star. In a follow-up study 
Ijerzvkiewicz & Sterken (1992) collected about one-hundred 
uvby observations in 1986 and found two significant peaks 
in the periodogram: at P — 3.02 d and at an alias pe- 
riod of 1.49 d. They also anal ysed the previous dataset of 



Vander Linden & Sterkenl(ll987l) . consisting of 22 observations 



between 1983 and 1984, and found the same periodicity, al- 
lowing them to improve the accuracy of the main period. 



Jerzvkiewicz & Sterkenl d 19921) concluded that the star may be 



an ellipsoidal variable, but that the 3.02 d period could not be 
ruled out to be caused by rotation or by high overtone g-mode 
pulsations. They found that rotational modulation is a reason- 
able working hypothesis for the star. So far there exists no pul- 
sation study or abundance analysis. The star has been found to 
be a photometric variable in the Hipparcos data (lKoen & Evei 
20021) with a frequency of 10.93482 d" 1 and amplitude 0.0031 



magnitude in V, based on a sinusoid fit to the 143 datapoints. 
The nature of this periodicity is not clear. 

In the course of our investigation of the role of magnetic 
fields in massive stars a number of indirect magnetic indicators 
have been identified dHenrichs et al. 20051) . The most signifi- 
cant indicator is a strictly periodic variation in the well-known 
UV stellar wind resonance lines, in particular C iv, Si iv, N v 
and Al m. The variation observed in these lines is, in these mag- 
netic objects, restricted to a velocity interval symmetric around 
the rest wavelength, in contrast to the often observed high- 
velo city absorption variations seen in m any O and B stars (see 



e.g. iKaper et al.l 1997; Fullertonl 120031) . This symmetric low 



velocity variability in the UV wind lines of B stars is encoun- 
tered among all studied magnetic oblique rotators (see exam- 
ples in Fig. [TJ which include many helium peculiar stars with 
spectral type around B1/B2 for the He-strong stars and around 
B7 for the He-weak stars. Prime examples are the He-strong 
B2V star HD 184927 where the var y ing U V lines (repeated in 
Fig. [TJ were discovered by Barker (1986), and the He weak 
B7 IIIp star a Scl, of w hich the rotation period was discovered 
by IShore et al.l dl990h from UV data (repeated in Fig. [TJ. To 
the group of early B-type stars with this same specific UV vari- 
ations also belong Cep (B1IV), V2052 Oph (B2IV-V) and 
( Cas (B2IV): respectively a /3 Cephei variable, a He-strong 
P Cephei star and an SPB (slowly pulsating B) star. The UV 
wind lines of these three stars and of the magnetic B stars 
were recognised to behave very similarly, which prompted a 
search for a magnetic field in these stars. After the prediction by 
Henrichs et al.l ( 1993 ), all three were indeed found to be mag- 



netic ( Henrichs et al 



2000ll2012tlNeiner et alJ2003albh . based 



on o bservations obtain ed with the Musicos spectropolarime- 
ter ( Donati et al.ll 19991) . The dipolar field strengths at the equa- 
tor of these stars are less than 1000 G. It should be noted that 
the very accurately determined rotation periods of these stars 
cover a rather wide range: 12.00 d, 5.37 d and 3.64 d for /? Cep, 
V2052 Oph and ( Cas, respectively. 



Gradv et al.1 dl987) already noted Civ variations in IUE 
spectra of cr Lup similar to those in f Cas (long before the 
magnetic origin of this behaviour was recognised), but it ap- 
pears now that the two spectra on which this conclusion was 



based were severely overexposed, inhibiting reliable quantita- 
tive results. 

No X-ray or radio flux from cr Lup has been detected, both 
indicators of a magnetic environment around the star. An upper 
limit of 10 3038 ergs -1 on the X-ray luminosity was obtained 



with the ROSAT satellite dBerghofer et alJll996l) . 

In this paper we report the analysis of the stellar wind lines 
of cr Lup and show that the variations observed in 4 subse- 
quently taken IUE spectra are without much doubt due to a 
corotating magnetosphere. In light of the above, we adopt the 
photometric period of 3.02 d as the rotation period of t he star , 
following the conclusion by Jerzvkiewicz & Sterkenl dl992 ). 
We reanalyse all available photometric data in Sect. [3] In Sect. 
4 we review the UV results. In Sects. [5]and[6]we report on mag- 
netic follow-up measurements of this star which confirmed our 
conjecture. In Sect. [7] we investigate the variations in equivalent 
width of selected spectral lines. The spectra allowed a subse- 
quent abundance analysis, which we present in Sect. [8] In the 
last section we summarise and discuss our conclusions. 

2. Stellar properties 



Acco rding to the Bright Star Catalogue ( Ho ffleit & Jasc hek 
199 lb cr Lup is a member of the local association Pleiades 
group. The adopted stellar parameters and other relevant ob- 
servational aspects are collected in Table [TJ We have adopted 
the re sults of the recent anal ysis based on high-quality spec- 
tra by Levenhagen & Leistei d2006l) . rath er than the ex trapo- 
lated and derived values as presented by iPrinia dl989h . who 
adopted a much lower effective temperature (17500 K) and 
therefore a much larger radius (8.8 R ) to be consistent with 
the lu minosity (which has a comparable value in the two pa- 
pers). ISokolovl d 19951) derived r eff = 25 980 ± 3250 K from 
the slope o f the Balmer continuuu m, consistent with the tem- 
perature of iLevenhagen & Leister (2006). They derived a pro- 
jected rotational veloc i ty of v sini = 80+14 kms -1 , whereas 
Brown & Verschuerenl ( 1997 ) found vsim = 69+10 kms 1 . 
Our analysis (Sect. [8j gives very similarly 68 ± 6 kms 1 , 



Table 1. Adopted and derived stellar parameters for cr Lup. 



Reference 



Spectral Type 


B1/B2V 


1 


V 


4.416 


2 


d(pc) 


176^ 


3 


log(L/L ) 


3.76 ± 0.06 


1 


T eS (K) 


23 000 + 550 


1 


logg (cms -2 ) 


4.02 ±0.10 


1 


R/Re 


4.8 ± 0.5 


Sect. [2] 


M/M e 


9.0 + 0.5 


1 


log Age (yr) 


7.13 ± 0.13 


1 


vsini (km s" 1 ) 


68 ±6 


Sect. [8] 


fphot (d) 


3.01938 ± 0.00022 


Sect. [3] 


Pmap, (d) 


3.01972 ± 0.00043 


Sect. [6] 


v rad (kms -1 ) 


0.0 ± 0.5 


Sect. H 



1. ILevenhagen & Leisteil d2006h - 2. SI MBAD - 3. Ivan Leeuwenl 
J2007l) - 4. Ijerzvkiewicz & Sterkerj < fl992h 
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Fig. 1. Gallery of magnetic B stars with signatures of magnetically-confined stellar winds in the C iv line. The typical modulation 
in the known magnetic B stars, the He-strong star HD 184927 (top left) and the He-weak star a Scl (bottom left) is very similar 
to that observed in the four stars /3 Cep, V2052 Oph, f Cas (see text for references) and cr Lup (this paper), which led to the 
discoveries of their magnetic fields. Note that rotation period, pulsation properties and helium peculiarity differ. In each figure 
the upper panel shows an overplot of all available International Ultraviolet Explorer (IUE) spectra, taken over many rotational 
cycles. The lower panel displays the ratio of the observed variation to the expected variation (due to the noise), showing the 
velocity range in the stellar rest frame within which significant variations occur. Note that the whole profile moves up and down, 
approximately symmetrically with respect to zero velocity. This phenomenon is only observed in magnetic early-type stars. 



which we adopt in this paper. With the stellar parameters in 
Table Q] we derive a radius of the star between 4.25 and 5.37 
R which is consistent w i th the radius of 4.5 R measured by 



Pasinetti Fracassini et al.l (12001b . This corresponds to the cal- 
culated value of P/sim of 3.0 ± 0.9 d, compatible with the main 
photometric period, which we identify therefore with the rota- 
tional period of the star. 

3. Photometric analysis and new observations 

Jerzvkiewicz & Sterkenl d 1992 ) determined the period P - 



3.0186 ± 0.0004 d from the two photometric datasets obtained 
as part of the Lo ng-term Photometry o f Variables (LTPV) pro- 
gramme at ESO dManfroid et aljl99ll) . If we define maximum 
light as phase zero and take the weighted average of all the re- 
ported epochs of maximum light in all four photometric colors 
of these two datasets, we find for the photometric ephemeris: 

7XLTPV) = HJD 2445603.65 ±0.11 

+ «x(3.0186 + 0.0004) (1) 



with n the number of cycles. Since iKoen & Evert (|2002) did 
not find this period in the Hipparcos data, we reconsid- 



ered the dataset. By using a CLEAN analysis ( Robe rts et al 



19871) we verified that the 3 d period is indeed not present. 
However, there are two obvious outliers at HJD 2448309.52 
and 2448818.90 which deviate more 86cr and 26<x, respec- 
tively. After removal of these two points, 158 measurements 



are left, spread over three years. A new CLEAN analysis re- 
vealed the strongest peak near 3.020 d, which is the same as 
the period from the ground-based observations. A least-squares 
fit gives PHippaicos = 3.01899 + 0.0006 4 d with amplitude 3.5 + 
0.4 mm ag, very similar to the values of Jerzvkiewicz & Sterken 
(1992). A combined fit of ground- and space-based pho- 
tometry, after normalising the average magnitude of the two 
datasets where we only used the Stromgren y (comparable to 
the Hp magnitude), gives as the best-fit photometric period 
^photometric = 3.01938 + 0.00022 d with amplitude 3.1 ± 0.5 
mmag. Maximum light occurs at: 



7\LTPV + Hipp) = HJD 2447620.48 + 0.56 
+ n x (3.01938 + 0.00022) 



(2) 



The best fit cosine function is overplotted on the phased 
data in Fig. [2] 

As an attempt to reproduce t he photometric light curve 
of IJerzvkiewicz & Sterkenl dl992l) . we conducted five nights 
of Stromgren v filter photometry of cr Lup, using the 0.9 m 
CTIO telescope operated by the SMARTS^ consortium with a 
2048 x 2048 CCD detector. Although cr Lup is a bright star, 
the associated errors, both instrumental and atmospheric, are 



1 The Small and Moderate Aperture Research Telescope System 
( http://www.astro.yale.edu/smarts/ 1 
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at or beyond the amplitude of the reported variability. We are 
thus unable with this instrumentation to contribute any further 
information regarding the photometric period. Furthermore, cr 
Lup lacks a suitable nearby comparison star, making differen- 
tial photometry a difficult task. 

4. UV observations 

An earlier UV observation with the OAO-2 satellite was re- 



HD 12738! a Lupi 
Wavelength (A) 



ported by IPanek & Sava ge (1976), who measured the equiva- 
lent widths of the Si iv and C iv doublets (4.5 A and 1.0 A). 

Table [2] contains a list of all available high-resolution UV 
spectra taken with the IUE satellite between 1983 and 1993. 
The first two spectra had to be discarded in our analysis, as 
they were severely overexposed, which prevented reliable re- 
sults in the regions of the stellar wind lines of interest because 
of the saturation and missing portions of the reduced spectra. 
Fig. [3] shows overplots of regions around the resonance lines 
of N v, Si iv, and C iv, along with the significance of their vari- 
ability. The C m complex near 1 175 A also shows similar vari- 
ability but in this wavelength region the signal-to-noise ratio 
(S/N) is too low to derive quantitative results. No other region 
in the spectrum between 1150 and 1900 A showed significant 
variability. We sampled the spectra upon a uniform grid with 
0.1 A spacing (degrading the resolution somewhat), and lo- 
cally transformed the wavelength into a velocity scale in the 
stellar rest frame, while correcting for the (small) radial veloc- 
ity of the star. The wavelengths of the strongest member of the 
doublet lines were taken as the zero point. We normalised the 
average flux of the whole spectrum by using only portions of 
the spectrum which were not affected by stellar wind variability 
or by echelle order overlap mismatches. The overall flux level 
did not differ more than 2% for the 4 spectra, likely within the 
expected photometric uncertainty of the instrument. We note 
that the absolute flu xes agree well with th ose measured with 
the OAO-2 satellite dCode & Meaddll979l) . The lower panels 
in the figures display the quantity cr b S /<x exp , which is the ra- 
tio of the observed to the expected value of sigma, the square 
root of the variance. The expected value of cr was calculated 
with a noise model as a fun ction of the flux / applicable for 
high-resolution IUE spectra (He nrichs et al Jl 19941) in the form 
of Atanh(//B) with best-fit parameters A — 19, representing 
the maximum signal to noise ratio, and B = 7.8 x 10~ 10 erg 

oLup B1/B2V 
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Fig. 2. LTPV with y filter (red dots) and Hipparcos (black 
squares) photometry as a function of phase with the best fit 
sinusoid overplotted. Estimated error bars are about 5 mmag. 
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Fig. 3. Overplot of the four selected profiles of N v (top), Si iv 
(middle), and Civ (bottom) of cr Lup. The two doublet rest 
wavelengths are indicated by vertical dashed lines. Top scale: 
wavelength. Bottom scale: velocity with respect to the stellar 
rest frame. In each panel the lower part displays the signifi- 
cance of the variability, expressed as the square root of the ra- 
tio of the measured to the expected variances. The changes are 
very similar to what is observed in known magnetic B stars (see 
Fig.Q]), but not in other B stars. 



cm _2 s _1 A _1 , representing the average of all flux levels outside 
the regions with variability (intrinsic or instrumental). 

All doublet lines show significant variations in the velocity 
interval [-200, +200] kms -1 . For all three ions the absorption 
line profiles in image number SWP 47374 are much deeper 
than in the other spectra, whereas the lines in SWP 47881 have 
the highest relative flux. This parallel behaviour excludes an 
instrumental effect as the origin. It exclusively occurs in mag- 
netic B stars, like the examples in Fig.Q] which made this star 
a prime candidate for magnetic measurements. 

Table [2] gives the equivalent-width measurements of the 
three wind lines over the velocity interval indicated. The listed 
values are the sum of the equivalent widths of the two doublet 
members. The continuum was taken at the averaged flux level 
of the two end points of t he interval. The error bars wer e cal- 
culated with the method of Chalabaev & Maillard| dl983l) . 
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Table 2. Journal of high resolution IUE observations of cr Lup, and the measured equivalent widths. 



Image 


Date 


HJD 1 


*exp 


Phase 2 


EW (NV) (A) 


EW (Si IV) (A) 


EW (C IV) (A) 


SWP 




-2440000 


S 




[-250, 250] kms" 1 


[-300, 300] kms" 1 


[-200,200] kms" 1 


20261 


1983-06-19 


5505.242 


109.6 


0.637 








22094 


1984-01-25 


5724.656 


99.8 


0.297 








45218 


1992-07-24 


8827.829 


45.6 


0.935 


0.163 + 0.073 


1.43 ± 0.11 


0.51 ±0.10 


47374 


1993-03-28 


9074.610 


49.8 


0.658 


0.469 ± 0.067 


1.96 ±0.10 


1.08 ± 0.09 


47881 


1993-06-16 


9155.283 


49.8 


0.374 


-0.060 ± 0.077 


1.32 ± 0.11 


0.16 ± 0.11 


48225 


1993-07-24 


9193.076 


47.8 


0.889 


0.095 ± 0.074 


1.46 ±0.11 


0.57 ± 0.10 



1. Calculated at mid exposure. - 2. Rotational phase, calculated with Eq.[5J 
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Fig. 4. Equivalent widths of the N v (top), Si iv (middle) and 
Civ (bottom) line profiles as a function of magnetic phase 
(Eq.0. 



The calculated phases of the IUE observations are given 
in column 5 of Table [2j where we used Eq. [5] below, in which 
phase corresponds to maximum (positive) magnetic field, and 
apparently also maximum light (see Sect. 16.21 ). From the figure 
it is obvious that all lines behave in the same manner, and that 
the two points near phase 0.9 have very similar values, which 
we interpret as a confirmation of the periodicity. Depending on 
the inclination angle i and the magnetic obliquity /3, a single 
or double sine wave with possibly unequal minima would be 
expected to fit the phased data, but with four points this is im- 
possible to verify. The lowest absorption (or highest emission, 
as the lines become filled in) appears around phase 0.4, i.e. 
compatible with minimum light. Based on the UV behaviour 
of other hot, magnetic stars, this would occur when a mag- 
netic pole points to the observer. This is opposite to what we 
would naively expect for a lightcurve with a single minimum. 
We expect that maximum absorption occurs when the magnetic 
equatorial plane, where most of the circumstellar material will 
be collected, intercepts the line-of-sight to the star. 



5. Magnetic observations and properties 

Initial spectropolarimetric observations in both left- and right- 
hand circularly polarised light were obtained at the Anglo- 
Australian Telescope (AAT) at six observing epochs, as shown 
in Table [3] Th ese were obtained using the SE MPOL spec- 
tropolarimeter JSemellll989T: iDonati et all 120031) on the AAT 
in conjunction with the University College London Echelle 
Spectrograph (UCLES). SEMPOL was positioned at the AAT's 
f/8 Cassegrain focus and consists of an aberration-free beam 
splitter and an achromatic A/4 plate orientated to transmit cir- 
cularly polarised light. SEMPOL splits the observed light into 
left-hand and right-hand circular polarisation and feeds it down 
twin fibres, which is fed into a Bowen-Walraven image slicer at 
the entrance to the high-resolution UCLES spectrograph. The 
A/4 plate can be rotated from +45° to -45° with respect to 
the beam splitter to alternate the polarisation in each fibre. In 
practice, four subexposures are obtained, alternating the wave 
plate orientation between +45°. These subexposures are com- 
bined in such a way as to eliminate spurious contributions to 
the polarisation spectrum. Further information on t he operation 



of the SEMPOL spec trop olarimeter is provided bvlSemel et al 



dl993h . lSemel dl989h and foonati et al 



(1997,2003). 



The detector used was a deep depletion EEV2 CCD with 
2048 x 2048 x 13.5/im square pixels. The UCLES was used 
with a 31.6 gr/mm grating covering 46 orders (#84 to #129). 
The central wavelength was set to ~526.8 nm with a coverage 
from ~437.7 nm to ~681.5 nm with a resolution of ~ 71000. 

Prompted by the magnetic detection in these spectra 
(see Fig. [5]), subsequent measurements were obtained with 
ESPaDOnS at the Canadian-France-Hawaii Telescope (CFHT), 
as part of the Magnetism in Massive Stars (MiMeS) Large 
Program. The ESPaDOnS spectropolarimeter and observing 
procedure are fundamentally similar to the SEMPOL/UCLES 
combination, providing circularly polarised spectra with a re- 
solving power R - 65000 from -370 to 1000 nm. A total 
of 6 AAT spectra and 20 CFHT spectra were acquired, with 
peak S/N pxL 1 ranging from 360-880 (AAT) and -500-1300 
(CFHT). For the log of observations see Table[3] 

Optimal extraction of the SEMPOL and ESPaDOnS echelle 
spectropolarimetric observations was done using the ESpRIT 
and Libre-ESpRIT (Eche lle Spectra Reduction: an Interactive 
Tool: IDonati et all 1997) reduction packages. Preliminary pro- 
cessing involved removing the bias and flat-fielding using a 
nightly master flat. Each spectrum was extracted and wave- 
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Table 3. Journal of observations and results of magnetic measurements of cr Lup. 



Nr. 


Obs. 1 


Date 2 


HJD J 


fexp 


S/N 4 


Phase 5 


B, b 


o-(B,) 


N,' 


aiNi) 








-2450000 


S 


pxr 1 




G 


G 


G 


G 


1 


A 


2007-12-28 


4463.24031 


4x400 


880 


0.141 


79 


20 


4 


20 


2 


A 


2008-01-01 


4467.25522 


4x400 


620 


0.470 


-69 


27 


-4 


27 


3 


A 


2008-12-15 


4816.25084 


4x400 


720 


0.043 


75 


25 


8 


25 


4 


A 


2008-12-16 


4817.24836 


4x400 


720 


0.373 


-39 


24 


-7 


24 


5 


A 


2008-12-17 


4818.24066 


4x400 


360 


0.702 


95 


49 





49 


6 


C 


2009-02-17 


4880.12491 


4x130 


1015 


0.195 


69 


24 


22 


24 


7 


C 


2009-02-17 


4880.13320 


4x130 


1188 


0.198 


96 


21 


-4 


21 


8 


c 


2009-02-17 


4880.14136 


4x130 


1140 


0.201 


52 


22 


9 


22 


9 


A 


2009-04-08 


4930.25755 


4x400 


750 


0.797 


57 


24 


-23 


24 


10 


c 


2010-02-27 


5255.10549 


4x200 


1114 


0.373 


-70 


22 


1 


22 


11 


c 


2010-03-02 


5258.09114 


4x200 


1021 


0.361 


-15 


24 


-37 


24 


12 


c 


2010-03-03 


5259.08881 


4x200 


779 


0.692 


-15 


31 


-15 


31 


13 


c 


2010-03-04 


5260.08395 


4x200 


1099 


0.021 


117 


22 


-20 


22 


14 


c 


2010-03-05 


5261.10469 


4x200 


959 


0.359 


-43 


25 


-6 


25 


15 


c 


2010-03-07 


5263.05040 


4x200 


1303 


0.004 


114 


19 


-13 


19 


16 


c 


2010-03-07 


5263.08299 


4x200 


1326 


0.014 


103 


19 


-3 


19 


17 


c 


2010-03-08 


5264.02698 


4x200 


490 


0.327 


-28 


50 


26 


50 


18 


c 


2010-03-08 


5264.03882 


4x200 


624 


0.331 


-30 


39 


27 


39 


19 


c 


2010-03-08 


5264.12734 


4x200 


744 


0.360 


-31 


33 


-40 


33 


20 


c 


2010-05-28 


5344.86081 


4x200 


908 


0.096 


73 


40 


3 


40 


21 


c 


2010-05-28 


5344.87315 


4x200 


579 


0.100 


14 


186 


28 


186 


22 


c 


2010-05-28 


5344.88541 


4x200 


664 


0.104 


173 


62 


75 


62 


23 


c 


2010-05-30 


5346.84258 


4x200 


980 


0.752 


27 


28 


23 


28 


24 


c 


2010-05-31 


5347.83115 


4x200 


1291 


0.079 


117 


22 


15 


22 


25 


c 


2010-06-01 


5348.88552 


4x200 


870 


0.428 


-123 


29 


52 


29 


26 


c 


2010-06-02 


5349.82460 


4x200 


963 


0.739 


1 


26 


11 


26 



1. A = obtained with the SEMPOL spectropolarimeter attached to the UCLES at the AAT; C = with ESPaDOnS attached to the CFHT. - 2. 
Date when the first observation started. - 3. Calculated at the center of the total exposure time (Y cxp ) of column 5. - 4. Quality of the Stokes V 
spectra, expressed as the signal to noise ratio per pixel around 550 nm in the reduced spectrum. - 5. Rotational phase, calculated with Eq.|5]- 
6. Longitudinal magnetic field values with their 1-cr uncertainties. - 7. Computed magnetic values of the diagnostic null (or N) spectrum with 
their l-<r uncertainties. 



length calibrated against a Th-Ar lamp, then normalised to the 
continuum using polynomial fits to individual spectral orders. 



5.1. Spectropolarimetric analysis using LSD 



As the Zeeman signatures in atomic lines are typically ex- 
tremely small (typical relative amplitudes of 0.1 per cent 
or less), we have ap plied the te chnique of Least-Squares 
Deconvolution (LSD. iDonati et al.lll997l) to the photospheric 
spectral lines in each echelle spectrum in order to create a sin- 
gle high-S/N profile for each observation. In order to correct for 
the minor instrumental shifts in wavelength space, each spec- 
trum was shifted to match the Stokes / LSD profi le of the tel- 
luric li nes c ontained in the s pectra , as described by Don ati et al 
(2003J) and iMarsden et al.l (|2006). The line mask used for the 
stellar LSD was a B2 star linelist containing 170 He and metal 
lines created from the Kurucz atomic database and ATLAS9 



atmospheric models (lKuruczlll993l) but with modified depths 
such that they are proportional to the depth of the observed 
lin es. Further informa tion on the LSD procedure i s provided 
bv lDonati et all d 1997b and lKochukhov et al.l d2010b . 



6. Spectropolarimetry: results 

The presence of a magnetic field on a star produces a signal in 
the Stokes V profile. The Stokes V profile is created by con- 
structively combining the left- and right-hand circularly po- 
larised light together from the 4 sub-exposures, while a diag- 
nostic null (TV) profile was created by destructively combin- 
ing the 4 sub-exposures. This null profile was used to check 
if any false detections have been introduced by variations in 
the observing conditions, instrument or star between the sub- 
exposures. Fig.[5]shows the Stokes V LSD profiles for the two 
discovery observations, which show opposite Zeeman signa- 
tures. As mentioned earlier, the Stokes V signature is much 
smaller than Stokes 7, typically less than 0.1% of the contin- 
uum. By combining the more than 170 lines, the S/N is im- 
proved dramatically. We used the same line list for all spectra. 
Table [3] includes the resulting mean SNR for each of the ob- 
servations. The observing conditions were often adverse, par- 
ticularly during the December 2008 run at the AAT as the star 
was very low in the East and was taken just prior to 12° twi- 
light. Hence the detections on the 15th and again on the 17th of 
December, 2008 were only marginal in LSD profiles sampled 
at the nominal resolution of the instrument. The observations at 
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o Lup B 1/B2V, magnetic data 28 Dec 2007 - 2 June 2010 
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Fig. 5. Discovery observations of the Zeeman signature in <x 
Lup at the AAT on December 28, 2007, (top) and four days 
later (bottom), at magnetic phases 0.14 and 0.47, respectively. 
Each figure shows the LSD Stokes unpolarised / (lower panel) 
and circularly polarised V (top panel) profiles of SEMPOL 
spectra, all normalised by the continuum. The middle panel dis- 
plays the (null) N profile for integrity purposes. The integrated 
Zeeman signature in the V profile over the width of the line 
(between the outer vertical dashed lines) gives the value of the 
longitudinal component of the magnetic field, integrated over 
the stellar surface. On the first night we obtained B[ = 79 + 20 
G, whereas on the second night the V signature changed sign, 
and we obtained B/ = -69 + 27 G. Note the irregular bumps in 
the / profiles, which indicate the presence of spots on the stellar 
surface. 



CFHT were all taken near airmass 3 because of the low decli- 
nation (-50°) of the star. The observation at JD 2455344.87415 
(number 21 in TableO was taken under very poor weather con- 
ditions, leading to a veer large uncertainties (null spectrum with 
N = 28 ± 186 G, and B = 14 ± 186 G). This observation was 
disregarded in further analysis. 



200 =- AAT + SEMPOL (red) 
150 r 
100 I 

50 ~- 
r 

- 5<,: \ 
-100 - 



It 



CFHT + Espadons (blue) 



-150 
4400 



4600 



4800 5000 
HJD - 2450000 



5200 



5400 



Fig. 6. Magnetic data as a function of time. Observations from 
AAT (red dots) and CFHT (blue squares) are indicated. 

6.1. Longitudinal magnetic field 

We derived the longitudinal component of the magnetic field, 
integrated over the visible hemisphere, from the Stokes / and V 
LSD profiles in the usual manner by calculating dMathvsll 989: 
Donaitiet alll997h : 



B f = (-2.14X 10 n G) 



j vV(v)dv 
Age f[l - I(v)]dv 



(3) 



where A, in nm, is the mean, S/N-weighted wavelength, c is the 
velocity of light in the same units as the velocity v, and g is 
the mean, S/N-weighted value of the Lande factors of all lines 
used to construct the LSD profile. We used A — 5 15 nm and g = 
1.27. The noise in the LSD spectra was measured and is given 
in Table [3] along with the S/N obtained in the raw data. The 
integration limits were chosen at + 90 kms~'. The longitudi- 
nal field is measured with a typical lcr uncertainty of 20-30 G, 
and varies between about +100 and -100 G. In contrast, the 
same quantity measured from the profiles show no signifi- 
cant magnetic field. The longitudinal field data as a function of 
time are shown in Fig. [6] and are reported in Table 3. 

6.2. Magnetic period and phase analysis 

A best fit to the magnetic data in Table [3]of the form 

Bf(t) = B + Bmax COS(27T(f /P + $)) (4) 

with l/cr(B) 2 error bars as weights yields P = 3.01972 + 
0.00043 d, B = 22 + 6G, B max = 96 + 8 G, and <f> = 0. 1 1 + 0.24 
with a reduced^ 2 = 1.01. An overplot with the data, folded in 
phase, along with the residuals is presented in the two lower 
panels of Fig. [7] The maximum magnetic field value occurs at: 



r(B max ) = HJD 2455 1 06 .0 1 ± 0.72 
+ n x (3.01972 + 0.00043). 



(5) 



The obtained period is very close to the photometric period 
of 3.01938 + 0.00022 d (see Eq.©. Extrapolating backwards to 
the epoch of photometric observations (with 2479 cycles), we 
find that the epoch of maximum light at HJD 2447620.48+0.56 
coincides within the errors with the epoch of expected maxi- 
mum magnetic field at HJD 2447620.86 + 0.72. As we have no 
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a priori knowledge of why maximum light would occur at max- 
imum field strength, we cannot derive a more accurate value of 
the rotational period from combining the phase information. 
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-0.5 
1.2 

1.1 



* 1 

a) 



0.9 
1.05 



aLup B1/B2V, Dec 2007- Jim 2010 

' 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 
L IUE data, 1992- 1993 , ' , 

f 



0.95 
1.2 



1.1 



* 1 

a) 



0.9 
200 
100 


-100 



U -200 
' 100 

5 50 

6 o 

"s -50 



aa 



? -100 



+ \ 



Al III 569i 




AAT + SEMPOL (red) 
CFHT + Espadons(blue) 

+ H t, 




Residuals 



-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 
Magnetic Phase (3.01972 d) 



Fig. 7. Magnetic and equivalent- width data (red dots: AAT; 
blue squares: CFHT) phased with Eq.[5] repeated over two ro- 
tational periods. Lower two panels: magnetic data with their 
residuals to the best-fit cosine curve. Panels 3, 4 and 5 from 
bottom: equivalent widths of the Si in 4575, He 1 4471 and Al m 
5696 lines with a best-fit sinusoid superposed. Note the oppo- 
site behavior of the different lines. Top: Civ 1540 equivalent 
width fifteen years earlier. The dashed curve is a suggested dou- 
ble sine wave, with arbitrary amplitude, and phased according 
to what is observed in other similar magnetic stars. 



6.3. Magnetic geometry 

The adopted values of the rotation period, vsim and stellar ra- 
dius (see Table [TJ imply that the inclination, i, of the rotation 



axis with the line of sight is about 58°, but not well constrained 
between 43° and 90°. For a centered dipole model, the mag- 
netic tilt angle with respect to the rotation axis, B, is then con- 
strained by the observed ratio B max /B m i n = cos(/ + B)/ cos(« - 
B) = -1.59+o^, implying B close to 90°. In other words, the 
line of sight is not so far (< 45°) from the rotational equato- 
rial plane, and per rotation we see both magnetic poles passing. 
Such a dipole geometry implies a polar field of B p 500 G. 



7. Equivalent widths 

The equivalent widths were measured for the following 15 
spectral lines: Cn ^4267, Hei A4471, 4920, 6678, Nn ,14601 
(blended with On), 4630, Onvl 4416, 4591, Sin ,15040, 5055, 
Sim ,14552 (blended with Nn), 4568, 4575, Aim 5696, and 
Ha /16563. The variation as a function of rotational phase of 
three representative lines are shown in Fig. [7] The three he- 
lium lines show the most significant modulation, in antiphase 
with the magnetic field, whereas opposite phase behavior, i.e. 
in phase with the magnetic field, is found among the five sil- 
icon lines. The slightly variable Al m varies in phase with the 
helium lines. The nitrogen lines were not detected to vary. None 
of the other studied lines show significant periodicity. This in- 
cludes Ha, which did not show any significant emission, see 
Fig. [8] We note that the antiphase behaviour between the He 
and Si lines is similar to what is observed in other magnetic 
He-strong and He-weak B stars, like for instance HD 37776, 
for which Khokhlova et al. (2000) showed that helium is over- 
abundant in regions where silicon is underabundant, and vice 
versa. The surface abundances of some of these elements are 
apparently not homogeneously distributed over the stellar sur- 
face, causing a rotational modulation of the equivalent widths, 
and in total luminosity. The irregular bumps in the / profiles in 
Fig. [5] are therefore likely explained by the presence of spots 
for a number of elements, associated with the magnetic field 
configuration. Zeeman Doppler Imaging techniques may give 
a more quantitative result. 



o Lup, Ha, CFHT 




6540 



6550 
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Fig. 8. Overplot of spectra around Ha, showing no obvious 
emission, and no significant variability. 
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8. Abundance determination 

For its spectral type it is not unusual that this star is magnetic: 
in fact it could be a He peculiar star. The ultraviolet He n line at 
1640.41 A is not prominent. Many magnetic early B stars have 
abundances deviating from solar. 

The available spectra include lines of H, He, C, N, O, Ne, 
Mg, Si, S, and Fe. The spectral lines of most of the elements 
are mildly variable, particularly in the shape and depth of the 
cores, but very little in the line wings. This strongly suggests 
that the atmospheric chemistry of <x Lup is not uniform, but is 
somewhat patchy. The fact that the overall equivalent width of 
most lines remains approximately constant in spite of the small 
core variations suggests that this patchiness is rather small scale 
- there is no strong evidence for variations from one magnetic 
hemisphere to the other, except for a mild variation of helium. 

To obtain a first characterisation of the atmospheric chem- 
istry of the star, mean abundances have been determined for 
several elements from each of two spectra. Such abundances, 
for elements that are certainly non-uniform on a small scale, 
correspond roughly to hemispheric average abundances for 
the hemisphere producing the observed spectrum. Two spectra 
have been analysed, one from near the negative field extremum 
and He line strength maximum at phase 0.428 (obtained on 
2010-06-02), and one from near field maximum and He mini- 
mum at phase 0.195 (obtained on 2009-02-17). 

The abundance determinations have been made using 
the line synthes is program ZEEMAN dLandstreet 1988; 
Wade et al.ll200ll) . designed to compute line profiles for mag- 
netic stars from a simple model of magnetic field structure 
(a low order multipole expansion) and a simple distribution 
model which in this case is taken to be a uniform distribu- 
tion over the star. It takes as input a model atmosphere in lo- 
cal thermodynamic equilibrium (LTE) computed with Kurucz' 
ATLAS for T e = 23000 K and logg = 4.0, and computes 
the observed spectrum. The program determines a best value 
of radial velocity v r = 0.0 + 0.5 kms , rotational velocity 
vsinz = 68+6 kms -1 , and abundance of a single element at 
a time b y optimising the fit to a segment of observed spectrum 
(see e.g. Bagn ulo et al ] l2004h . iterating this process until a sat- 
isfactory fit. A dipolar field with polar strength = 400 G 
observed at about 45° has been assumed, but the results are 
essentially independent of the assumed field geometry. 

The program assumes LTE. Although this is not a good 
approximation in ge neral for a star of T e - 23000 K, 
Przvbilla et al.l d201 ll) showed that LTE modelling can yield 
meaningful abundances up to about this effective temperature, 
provided the right spectroscopic indicators are employed. They 
provide a list of some of the indicators to use. We have iden- 
tified others by synthesising an available spectrum of a Pyx 
= HD 74575, a star very similar to cr Lup, for which the 
fundamental pa rameters and abunda nce table are known with 
high accuracy (IPrzvbilla et al.l 120081) . By creating a synthetic 
spectrum of a Pyx based on these stellar parameters, we can 
identify spectral lines whose strengths are correctly computed 
in LTE, and then use these lines to determine abundances in 
cr Lup. 



Table 4. Characteristic chemical abundances of <x Lup, de- 
termined at rotational phases near magnetic field extrema 
(columns 2 and 3), and compared with solar values (column 
4). 

Element log(^) log(^) log(^) [X/H] 

°\"H/0.195 c V"II/0.428 °\"ll/o 



He 


-0.88 ±0.12 


-0.75 


±0.14 


-1.07 


+0.19 to +0.32 


C 


-3.90 + 0.19 


-3.85 


±0.15 


-3.57 


-0.30 


N 


-3.73+0.16 


-3.75 


±0.15 


-4.17 


+0.43 


O 


-3.28 + 0.16 


-3.33 


±0.14 


-3.31 


+0.00 


Mg 


-4.57 + 0.15 


-4.52 


±0.15 


-4.40 


-0.15 


Al 


-5.80 + 0.15 


-5.75 


±0.15 


-5.55 


-0.22 


Si 


-5.0 + 0.4 


-5.0 


±0.4 


-4.49 


-0.5 


S 


-5.68 + 0.15 


-5.68 


±0.15 


-4.88 


-0.80 


Fe 


-4.60 + 0.10 


-4.60 


±0.10 


-4.50 


-0.10 



Abundances were determined using several ~100 A win- 
dows. A sample of the fit for phase 0.428 is shown in Fig. [9] 
The results are summarised in Table|H which lists for each ele- 
ment studied both the logarithmic value of the ratio of the num- 
ber density of the element to the number density of H for each 
of the two spectra , the v alue of this ratio in the Sun as reported 
bv lAsplund et al. (2009), and the logarithmic ratios of the val- 
ues in cr Lup relative to the Sun (average values for elements 
that have closely similar abundances in the two spectra). 

This table needs several remarks. Uncertainties in abun- 
dance are estimated from the scatter in the value derived from 
different spectral windows. This scatter probably arises in part 
from inaccuracy in atomic data, in part from probably not pre- 
cisely correct fundamental parameters, in part from neglect of 
non-LTE effects, and in part from the small-scale abundance 
inhomogeneities that are not modelled. 

For He, we have used the relatively weak He i lines at 4437, 
4713, 5015, and 5047 A, which do not have signficant forbid- 
den components, and for which LTE an d non-LTE calculation s 
give very similar profiles at 22000 K JPrzvbilla et al.l l20lll) . 
We also modelled the very weak He i line at 4686 A; this line 
is approximately the correct strength when computed with the 
abundance of He derived from the He i lines, consistent with 
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Fig. 9. Portion of the spectrum at phase 0.428, with most lines 
identified, overplotted with the model fit to determine abun- 
dances. 
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the values of T e and log g in Table Q] The He abundance in 
both spectra is weakly overabundant with respect to the solar 
value, by about a factor of two, so cr Lup is a mild He-rich star. 
An abundance difference of about 0.1-0.15 dex is consistently 
found between the two spectra for all the lines modelled, with 
a larger abundance in the spectrum at phase 0.428, consistent 
with the larger equ i valen t width observed for He 1 447 1 (Pig. |9j>. 
Przvbilla et al.1 ( 201 lh point out that many weak lines of N n 



and O n are well described by LTE, but that C n lines are mostly 
out of LTE. For C they recommend using the lines of C n mul- 
tiplet (16) around 5145 A for LTE abundance determination. 
Accordingly, the abundances of N and O are based on several 
weak lines of each element in several spectral window between 
4400 and 5050 A, while those of C are based on the recom- 
mended multiplet. The uncertainty is estimated from the dis- 
persion of the abundances found in various spectral windows. 
From these data we find that there is no significant difference 
in the abundance of these three elements between the two mag- 
netic hemispheres. C is found to be underabundant with respect 
to the Sun by about a factor of 2; N is overabundant by about a 
factor of 2, and O is approximately solar. Abundances of Mg, 
Al, S, and Fe are determined using weak lines of the dominant 
ionisation states. None seems to be significantly variable, and 
are all mildly subsolar in abundance, except for S which is al- 
most one dex underabundant. 

A number of lines of Si n and m, and a few lines of Si iv, 
are visible in the spectra of a Pyx and cr Lup. When these 
Si lines are synthesised for a Pyx in LTE using the parame- 
ters and abundances of this star, most of the observed lines are 
found to be very discordant with the computed spectrum. The 
same problem is present in the spectrum of cr Lup. The abun- 
dance of Si in cr Lup was initially obtained using the weak 
Si ii lines at 5040 and 5056 A, and the considerably stronger 
Si m triplet 4552, 4567, and 4574 A. The two levels of ioni- 
sation gave abundances that differed by about a factor of 30, 
from about -5.50 for the Si n lines to -3.95 for the Si m lines. 



From the discussion of this element bv lPrzvbilla et al.l (1201 lh . 
it appears that this is largely a problem of major departures of 
Si levels from LTE, an d that even non-LTE calculations (e.g. 
Becker & Butlerlfl990l) do not fully account for the departures. 

Accordingly, we have simply searched the spectrum of 
a Pyx for lines of Si t hat are approximately correct with the 
Si abundance found by iPrzvbilla et al.l (12008b . The most use- 
ful lines seem to be Si iv at 4088 and 41 16 A. These lines are 
consistent with an abundance of Si that is about a factor of 4 
below the solar value, in both the spectra modelled. However, 
we give this element an unusually large uncertainty in view of 
the exceptional difficulties of LTE analysis. Although the exact 
value is extremely uncertain, it appears that the abundance of 
Si in cr Lup is probably no higher than the solar ratio. As ex- 
pected, line profiles of most Si n and m lines computed with 
this abundance do not match the observed profiles at all well. 

9. Discussion and conclusions 

In this paper we predicted and confirmed that cr Lup is an 
oblique magnetic rotator. The period determined from photo- 
metric measurements is 3.01938 + 0.00022 d, while that deter- 



mined from the magnetic measurements is 3.01972+0.00043 d. 
These periods are in formal agreement, and we interpret them 
to be the rotational period of the star. Using the +100 G varia- 
tion of the longitudinal magnetic field, in combination with our 
estimate of the rotation axis inclination i 58°, we model the 
magnetic field as a dipole with obliquity [} - 90° and surface 
polar field strength Z?d — 500 G. Modeling of the spectrum us- 
ing the LTE synthesis code ZEEMAN reveals enhanced chem- 
ical abundances of a variety of elements. We observe variations 
of the equivalent widths and shapes of line profiles that repeat 
with the rotation period, indicating that chemical elements are 
distributed non-uniformly across the stellar surface, with pat- 
terns that depend on the element in question. 

These properties do not stand out in any particular way 
from the general population of early B-type magnetic stars. 
The longitudinal field variation is approximately sinusoidal, in- 
dicative of an important dipolar component to the field. The 
polar field strength is not atypical for such stars, nor is the 



large inferred value of the obliquity (see, e.g. Peti t et alJl201 1 



Alecian et al.l l201 1). The rotational period is within the range 
of periods observed for such stars, and the chemical peculiari- 
ties and their non-uniform distributions are similar to those of 
other stars with similar temperatures. We therefore conclude 
that cr Lup is a rather n ormal magnetic B-type star. 

As summarised bv lJerzvkiewicz & Sterkenl ( 11992b . similar 
phot ometric periods have been found in a number of Be and Bn 



stars dBalona et al.ll987 : van Vuuren et alJl988T:ICuvpers et al. 
1989b iBalona et al.l Il99lh . see also ICarrier & Burkil d2003l) . 
termed A Eri variables. These are short-term periodic Be stars 
with strict ly periodic var i ations with periods in the range 
0.5-2.0 d. Vogt & Penrod (1983) suggested that these vari- 
ations could be explained either by non- radial pulsation s or 
by clouds close to t he stellar photosphere. Balona| dl99ll) and 
Balo na et al.l dl991) suggested a rotational origin, specifically 
caused by material trapped above the photosphere by magnetic 
fields. I n this aspect it is interesting to mention that lNeiner et al] 
(2003c) reported a possible detection of a magnetic field in the 
Be star a> Ori, which also belongs to this class (alt hough recent 
result s are unable to confirm the repor ted field: see|Neiner etaD 
( 2012 ) . As was demonstrated by e.g. Krticka et al] "( 2009 ) and 
assumed in this paper, spots caused by surface abundance inho- 
mogeneities, associated with the magnetic field, cause the peri- 
odic photometric variations. Further investigation is needed to 
establish whether a similar phenomenon holds for A Eri stars as 
well. 

This star is the fourth confirmed B star which has been di- 
agnosed to be magnetic from the variability of its ultraviolet 
stellar wind lines, i.e. together with (3 Cep, V2052 Oph and ( 
Cas. These stars, however, had their rotation period established 
from UV data, whereas the period of cr Lup was found from 
photometry. It could be due to the short rotation period that 
the photometric changes are more easily detectable in this star, 
and/or that surface spots are less pronounced in the other stars. 
From these four B stars the presence of a magnetic field appar- 
ently does not depend on the pulsation or rotation properties. 

It is interesting to com pare the abundance s of these mag- 
netic B stars, as reported bv lMorel et al.l d2006l) . All these stars 
are enhanced in nitrogen, whereas all but /? Cep are helium 
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enriched. Apparently nitrogen enrichment in the presence of 
a magnetic field is even more common than helium over- or 
under-abundance. This property was in fact used to search for 
a magnetic field in ( Cas. In this context we also note that in 
the previously-known nitroge n enriched B0.5V s tar £ 1 CMa a 
magnetic fiel d was found by [ Hubrig et alJ ( 2006h . which was 
confirmed bv ISilvester et alJ (120090 . On the other hand, other 
some nitrogen e nriched early B stars (e.g. 5 Cet) are not found 
to be magnetic dWade & Mimes Collaboration 201 1 ). Thus N 
enrichment, while not necessarily a direct tracer of magnetism 
in B stars, could still serve as an additional indicator that will 
be useful to discover more magnetic B stars. 
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